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HIGHLIGHTS 


•  The  effects  of  CO  on  spatial  PEMFC  performance  and  EIS  were  studied  with  a  segmented  cell. 

•  The  experiments  were  performed  using  various  cathode  gases:  air,  O2  and  H2. 

•  Injection  of  CO  resulted  in  a  voltage  decrease  and  redistribution  of  segments'  currents. 

•  The  spatial  EIS  data  were  analyzed  using  the  equivalent  electric  circuits  approach. 

•  A  current  distribution  model  and  EIS  interpolation  method  were  applied  for  detailed  analysis. 
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This  paper  presents  experimental  and  modeling  results  of  the  effect  of  low  CO  concentration  (2  ppm)  on 
the  spatial  performance  of  PEMFC  as  well  as  its  spatial  electrochemical  impedance  spectroscopy  (EIS) 
responses.  The  cell  was  operated  at  constant  current  using  various  cathode  gases:  air,  02  and  H2.  Due  to 
CO  adsorption  on  the  Pt  anode  and  its  poisoning,  the  cell  voltage  decreased  and  spatial  current  redis¬ 
tribution  was  observed.  The  steady  state  voltage  losses  were  0.089,  0.280  and  0.295  V  for  the  H2/O2,  H2/ 
air  and  H2/H2  gas  configurations,  respectively.  EIS  data  revealed  a  pseudo-inductive  behavior  in  the  low 
frequency  region  for  inlet  segments  of  the  cell  operated  under  H2/air  and  H2/H2  conditions.  Operation 
with  O2  as  an  oxidant  did  not  cause  any  pseudo-inductance.  Analysis  of  the  EIS  and  anode  overpotential 
data  suggested  that  CO  oxidation  occurred  via  chemical  or  electrochemical  mechanisms,  or  a  combi¬ 
nation  of  both  depending  on  the  selected  cathode  gas.  The  spatial  EIS  data  were  analyzed  using  the 
equivalent  electric  circuits  approach.  The  distributions  of  the  equivalent  electric  circuit  parameters  are 
presented  and  discussed.  A  current  distribution  model  and  EIS  interpolation  technique  were  successfully 
applied  for  detailed  analysis  of  CO  effects  on  the  spatial  PEMFC  performance  and  EIS. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cells  are  receiving  increasing  attention  due  to  their  high 
energy  conversion  rates  and  harmless  emission  products  as  well  as 
their  promising  future  applications  for  powering  stationary  and 
portable  devices  and  electric  vehicles  [1  ].  The  highest  performance 
is  achieved  with  H2,  which  is  the  preferred  fuel  for  low  temperature 
proton  exchange  membrane  fuel  cells  (PEMFCs).  However,  the 
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application  of  H2  as  a  fuel  for  PEMFCs  has  several  limitations, 
predominantly  due  to  current  methods  of  H2  production  which 
commonly  occurs  via  the  stream  reforming  of  natural  gas  [2,3]  as 
well  as  water  electrolysis.  Focusing  on  the  reforming  processes,  a 
by-product  of  these  is  carbon  monoxide.  To  minimize  the  CO  con¬ 
tent,  several  additional  stages  are  incorporated  into  the  reforming 
process,  namely  the  water  gas  shift  reaction  and  low  temperature 
preferential  oxidation.  The  resulting  CO  concentration  is  in  the 
range  of  10—50  ppm.  However,  even  a  low  CO  content  can  signifi¬ 
cantly  reduce  the  performance  of  fuel  cells  [4—17],  Recent  stan¬ 
dards  released  specifying  acceptable  hydrogen  gas  compositions 
for  fuel  cell  applications  set  the  maximum  CO  concentration  at 
0.2  ppm  in  attempt  to  alleviate  concerns  of  contamination  (ISO 
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14687-2,  SAE  J2719).  The  main  reason  for  the  performance  loss  due 
to  CO  contamination  in  PEMFCs  is  the  preferential  adsorption  of  CO 
on  Pt  sites,  which  leads  to  the  inhibition  of  H2  dissociative 
adsorption  and  its  further  oxidation.  Several  approaches  have  been 
proposed  to  mitigate  the  effect  of  CO  on  PEMFCs.  These  approaches 
include  preventative  measures  by  reducing  carbon  monoxide 
concentrations  through  current  advances  in  preferential  CO 
oxidation  techniques  and  the  water  gas  shift  reaction  [1—3,18],  and 
in-situ  techniques  such  as  the  pulsed  oxidation  of  a  CO  poisoned 
electrode  [19,20],  air  bleeding  into  a  H2  stream  [3,4,13,21—25],  the 
use  of  CO-tolerant  catalysts  [7,10,11,23,26,27],  and/or  operation 
under  particular  conditions  (high  temperature,  pure  H2  pulsing, 
high  cathode  back  pressure)  [28—31], 

For  the  diagnosis  of  fuel  cell  performance,  electrochemical 
impedance  spectroscopy  (EIS)  is  a  well-established  method 
[32—35],  The  method  is  based  on  a  small  sinusoidal  perturbation  of 
the  electrochemical  system  (either  in  current  or  voltage)  and 
measurement  of  the  response  signal.  Because  of  the  small  ampli¬ 
tude  of  the  perturbation,  the  system  is  assumed  to  be  linear  in  the 
range  of  perturbation,  and  the  linear  impedance  is  analyzed.  Vari¬ 
ation  of  the  perturbation  frequency  allows  processes  with  different 
time  constants  to  be  separated.  EIS  has  been  applied  for  the 
investigation  of  various  phenomena  in  PEMFC:  dehydration, 
flooding  and  poisoning  by  impurities  including  CO. 

The  first  in-situ  measurements  of  CO  poisoning  with  EIS  were 
published  by  B.  Muller  [36]  and  M.  Ciureanu  [37],  The  paper  [36] 
reported  the  time  dependence  of  impedance  spectra  under  gal- 
vanostatic  conditions  and  CO  exposure  of  the  Pt  anode,  while 
potentiostatic  electro-oxidation  of  H2/CO  was  studied  in  Ref.  37], 
It  is  useful  to  operate  the  fuel  cell  in  the  galvanostatic  mode 
because  a  constant  current  density  means  a  constant  reaction  rate. 
Additionally,  the  impedance  of  the  membrane  and  cathode  might 
also  be  assumed  to  be  constant,  and  any  observed  differences  in 
EIS  profiles  can  be  attributed  to  changes  in  the  anode  impedance. 
Several  other  publications  also  present  detailed  EIS  data  con¬ 
cerning  the  effect  of  CO  on  H2  oxidation  at  the  anode  [38—50],  It 
was  demonstrated  that  the  CO  poisoning  of  Pt-containing  anodes 
caused  an  increase  of  the  total  impedance  of  a  cell.  The  impedance 
response  strongly  depended  on  the  exposure  time,  potential,  CO 
concentration,  temperature,  and  oxygen  crossover.  An  occurrence 
of  a  pseudo-inductive  loop  was  detected  at  low  frequency  under 
certain  conditions.  The  loop  was  observed  to  be  related  to  the 
surface  relaxation  of  the  anode.  Furthermore,  the  effect  of  CO  on 
the  electrochemical  behavior  of  the  fuel  cell  can  be  interpreted 
using  the  Faraday  impedance  in  addition  to  the  potential- 
dependent  hindrance  of  the  charge  transfer 
[36,38—40,43,46,51—53],  Progressive  poisoning  with  CO  of  a  fuel 
cell  during  EIS  measurements  can  complicate  the  analysis  of  the 
data  because  the  examined  system  should  be  at  steady-state  when 
EIS  is  running.  However,  violation  of  steady-state  conditions  often 
occurred,  and  to  solve  these  problems,  time-resolved  EIS  (TREIS) 
was  introduced  [54], 

The  use  of  H2  fuel  containing  low  levels  of  CO  is  likely  to  cause  a 
varying  distribution  of  Pt  sites  occupied  by  CO,  thus  creating  areas 
with  low  performance  located  mainly  at  the  inlet  of  the  cell, 
leading  to  the  redistribution  of  local/spatial  currents  and  affecting 
the  spatial  impedance  responses.  Therefore,  EIS  of  the  cell  under 
CO  exposure  depends  not  only  on  the  time  of  exposure  but  also  on 
the  location  inside  the  membrane  electrode  assembly  (MEA).  The 
evaluation  of  fuel  cell  performance  with  a  single  cell  approach  only 
provides  an  average  of  the  local  voltage,  current,  and  impedance 
values  and  does  not  reveal  the  spatial  behavior  of  the  cell.  A 
segmented  cell  system  is  a  powerful  tool  for  understanding  the 
details  of  locally  resolved  fuel  cell  processes  [55],  There  are  many 
examples  of  segmented  cell  research  applications  in  PEMFC 


studies,  including  basic  investigations  of  local  current  distributions 
[56— 60  ,  gas  and  water  management  effects  [61—65  ,  stack  and 
single  cell  diagnostic  techniques  [66—69],  defect  detection  and 
localization  methods  [70—74],  recirculation  [75],  start-up,  and 
starvation  impact  [76—79  .  However,  there  are  only  a  few  papers 
devoted  to  using  segmented  cells  with  poisoned  electrodes 
[16,17,22,30,80],  The  authors  of  these  publications  discussed  the 
redistribution  of  the  local  currents  under  CO  poisoning  at  rela¬ 
tively  high  concentrations  on  the  order  of  50  ppm  up  to  3%.  Our 
previous  paper  [6  presented  results  concerning  the  effect  of  low 
concentration  CO  (2  ppm)  on  PEMFC  spatial  performance  under 
galvanostatic  conditions  (0.8  A  cm~2)  and  H2/H2,  H2/air  and  H2/O2 
gas  configurations  for  the  anode/cathode.  In  that  paper,  spatial  EIS 
data  presented  the  dynamic  responses  of  the  segments  upon  CO 
exposure  and  revealed  pseudo-inductance  at  low  frequency  and  an 
increase  not  only  in  the  anode  charge  transfer  resistance  but  also 
in  the  cathode  charge  transfer  resistance.  This  paper  is  a  contin¬ 
uation  of  our  studies  and  focuses  on  detailed  spatial  EIS  charac¬ 
terization  of  PEMFCs  during  low  CO  concentration  exposure 
(2  ppm).  The  selected  CO  concentration  is  still  1  order  of  magni¬ 
tude  higher  than  current  H2  fuel  cell  commodity  gas  specifications, 
but  in  comparison  to  other  studies  mentioned  is  quite  low  and  for 
the  purposes  of  this  analysis  may  be  considered  to  more  closely 
represent  the  real  operating  conditions  for  PEMFCs.  Moreover,  a 
mathematical  model  describing  the  temporal  evolution  of  the 
spatial  current  densities  was  adapted  and  expanded;  results  are 
presented  and  discussed  to  help  further  interpretation  of  the 
experimental  observations. 

2.  Experimental 

All  the  experiments  were  conducted  on  a  single  cell  test  station 
using  Hawaii  Natural  Energy  Institute’s  (HNEI)  segmented  cell 
system,  which  enables  the  simultaneous  acquisition  of  spatially 
distributed  data  [68],  The  segmented  cell  approach  for  this  study 
builds  upon  the  works  of  Cleghorn  et  al.  [56],  the  German  Aero¬ 
space  Centre,  Stuttgart  [81],  Ballard  Power  Systems  Inc.  [57],  and 
Los  Alamos  National  Laboratory  (LANL)  [82],  HNEI's  segmented  cell 
system  is  partially  based  on  the  LANL  design  using  closed  loop  Hall 
sensors  and  an  improved  data  acquisition  system.  These  en¬ 
hancements  allow  the  system  to  perform  simultaneous  rather  than 
sequential  measurements  of  spatial  EIS,  spatial  linear  sweep  vol¬ 
tammetry  (LSV),  and  cyclic  voltammetry  (CV).  The  segmented  cell 
hardware  is  based  on  the  HNEI  100  cm2  cell  design.  The  hardware 
contains  a  segmented  anode  flow  field  consisting  of  ten  consecu¬ 
tive  segments  disposed  along  the  path  of  a  ten-channel  serpentine 
flow  field.  Each  segment  has  an  area  of  7.6  cm2,  with  its  own 
distinct  current  collector.  The  same  channel  designs  are  used  for 
both  the  segmented  anode  and  the  standard  cathode  flow  fields 
(the  reactant  streams  were  arranged  in  a  co-flow  configuration). 

The  segmented  cell  system  consists  of  the  cell  hardware,  the 
current  transducer  system,  the  data  acquisition  device  and  a  single 
cell  test  station.  The  current  transducer  system  was  custom 
designed.  A  closed  loop  Hall  sensor  (Honeywell  CSNN  191)  is 
employed  for  current  sensing.  For  EIS  measurements  these  sensors 
show  very  little  inductance  over  the  entire  frequency  range  of  in¬ 
terest.  The  system  allows  the  investigation  of  as  many  as  10  current 
channels  in  a  high  (standard)  current  mode  and  16  channels  in  a 
low  current  mode.  The  standard  current  mode  enables  the  mea¬ 
surement  of  segment  currents  up  to  15  A.  The  current  limit  of  the 
data  acquisition  system  can  be  extended  to  30  A  or  more  using  a 
unique  counter  current  technology  that  allows  a  flexible  segmented 
cell  design,  high  current  operation  and  increased  accuracy  during 
EIS  experiments.  The  low  current  mode  of  the  system  yields  current 
measurement  up  to  375  mA,  with  an  accuracy  of  ±2.5%,  which  is  a 
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convenient  feature  for  CV  or  LSV  experiments.  Voltage  and  current 
signal  data  collection  was  performed  with  a  National  Instrument 
PXI  data  acquisition  instrument  operating  on  HNEI-developed 
LabView  programs. 

The  cell  hardware  was  operated  with  standard  100  cm2  Ion 
Power  MEAs.  The  anode  and  cathode  electrodes  were  composed  of 
a  50%  Pt/C  catalyst  coated  on  a  DuPont  Nation  212  membrane  with  a 
loading  of  0.4  mgPt  citT2.  We  used  a  Pt/C  catalyst  as  the  model 
system  in  this  work  because  it  exhibits  a  noticeable  response  even 
at  low  CO  poisoning  concentrations.  The  gasket  material  was  made 
of  Teflon,  with  thicknesses  of  203  pm  and  178  gm  for  the  anode  and 
cathode,  respectively.  Sigracet  25  BC  was  used  as  the  anode  and 
cathode  GDLs.  The  anode  used  a  segmented  GDL  and  gasket 
configuration,  whereas  a  single  GDL  was  applied  at  the  cathode  and 
an  active  area  of  MEA  was  76  cm2. 

During  CO  exposure,  the  dry  contaminant  gas  was  injected  into 
the  humidified  feed  stream  of  the  anode,  keeping  the  humidifica¬ 
tion  of  the  gas  constant  by  increasing  the  temperature  setting  of  the 
humidifier  unit.  The  anode/cathode  testing  conditions  for  the 
contamination  experiments  were  H2/air  (or  O2,  or  H2),  48.3/ 
48.3  l<Pag  backpressure,  100/50%  relative  humidity  and  2/2  stoi¬ 
chiometry.  The  cell  temperature  was  60  °C.  The  cell  was  operated  at 
a  current  density  of  0.8  A  cm-2  (based  on  the  76  cm2  of  the  active 
MEA  area).  The  contaminant  experiment  was  divided  into  three 
phases: 

•  Phase  1  was  a  pre-poisoning  period  conducted  with  neat  H2.  The 
minimum  length  of  Phase  1  was  approximately  1/3  of  the  length 
of  the  entire  experiment. 

•  Phase  2  was  a  poisoning  period  conducted  with  a  constant 
amount  of  contaminant  injected  into  the  anode  feed  stream. 

•  Phase  3  was  a  recovery  period  conducted  again  with  neat  H2. 

The  contamination  experiment  was  conducted  in  H2/air,  H2/02 
and  in  a  driven  H2/H2  mode.  For  the  H2/H2  test,  the  gas  flow  rate 
was  0.847  1  min-1  and  the  relative  humidity  was  100%  for  both  the 
anode  and  cathode. 

The  polarization  curves  (VI  curves)  in  the  H2/air  gas  configura¬ 
tion  were  measured  under  the  same  conditions  as  the  contamina¬ 
tion  experiment.  VI  curve  measurements  were  also  performed  in 
the  H2/02  configuration.  In  the  latter  case,  the  same  flow  rates  were 
used  as  in  the  H2/air  configuration  to  keep  the  water  transport  in 
and  out  of  the  cell  identical  at  a  given  total  cell  current  density. 
Consequently,  the  stoichiometry  for  02  increased  to  9.5.  The 
resulting  two  different  VI  curves  (H2/air  and  H2/02)  were  used  for 
the  determination  of  a  segment's  activation,  ohmic  and  mass 
transfer  overpotentials,  as  described  previously  [68], 

The  VI  curve  measurements  and  contamination  experiments 
were  combined  with  EIS  to  determine  the  cell's  and  segments’ 
high-frequency  resistances  (HFR)  and  to  measure  the  electro¬ 
chemical  impedance  spectra  for  all  ten  segments  and  the  overall 
cell.  The  selected  frequency  range  for  the  EIS  experiments  was 
0.05  Hz— 10  kHz  (some  measurements  were  performed  in  the  fre¬ 
quency  range  of  0.01  Hz— 10  kHz),  and  the  amplitude  of  the  sinu¬ 
soidal  current  signal  perturbation  was  2  A,  which  resulted  in  a  cell 
voltage  response  of  10  mV  or  lower.  The  HFR  was  determined  from 
the  intercept  of  the  EIS  with  the  x-axis  at  higher  frequencies.  The 
resulting  impedance  spectra  were  analyzed  based  on  the  fitting  of 
an  equivalent  electrical  circuit  (EEC)  model.  The  values  of  the  circuit 
elements  were  determined  by  fitting  the  experimental  data  using 
ZView  software  (Scribner  Associates)  at  0.05  Hz  to  10  kHz. 

CV  experiments  were  conducted  for  the  determination  of  the 
electrochemically  active  surface  area  (ECA)  using  a  Parstat  2273 
potentiostat/galvanostat  from  Princeton  Applied  Research.  CVs 
were  performed  at  a  cell  temperature  of  35  °C  with  a  scan  rate  of 


20  mV  s  ,  while  100%  humidified  hydrogen  and  nitrogen  were 
supplied  to  the  reference/counter  and  working  electrodes,  respec¬ 
tively,  at  a  flow  rate  of  0.750  1  min-1.  For  each  measurement,  three 
cycles  were  applied  over  a  potential  range  from  -0.015—1.1  V  vs. 
the  hydrogen  reference  electrode  (HRE).  The  hydrogen  desorption 
peak  area  of  the  third  cycle  was  used  to  determine  the  ECA. 
Hydrogen  crossover  experiments  were  performed  at  the  same 
temperature  and  flow  conditions  as  the  ECA  experiments  using  a 
single  potential  sweep  from  0.1  to  0.4  V  vs.  the  HRE  at  a  scan  rate  of 
0.1  mV  s_1. 

Before  the  CO  exposure  experiment,  the  segmented  cell  was 
assembled  using  established  procedures,  conditioned,  and  sub¬ 
jected  to  a  set  of  beginning  diagnostic  tests.  The  diagnostics  con¬ 
sisted  of  CV  to  determine  the  ECAs  of  the  electrodes,  LSV  to 
determine  the  hydrogen  crossover  current  of  the  cell,  and 
measuring  of  VI  curves  with  spatial  EIS  using  the  H2/air  and  H2/02 
gas  configurations  to  determine  the  performance,  overpotential 
distributions  of  the  cell,  and  impedance  responses.  The  same  set  of 
diagnostic  tests  was  repeated  after  the  CO  experiment.  A  compar¬ 
ison  of  the  results  before  and  after  the  CO  exposure  provided 
insight  into  the  effects  of  the  contaminant  on  the  spatial  properties 
of  the  cell.  The  same  MEA  sample  was  used  for  three  CO  exposure 
tests  performed  at  different  gas  configurations:  H2/02,  H2/air  and 
H2/H2.  The  MEA  was  reconditioned  after  each  test. 


3.  Modeling  of  CO  poisoning 

3.1.  Equivalent  electric  circuits  models  for  analysis  of  obtained  EIS 
data 

EECs  which  have  been  found  suitable  for  fitting  EIS  data  without 
CO  poisoning  are  presented  in  Fig.  1.  The  EECs  have  several  main 
parts: 

1 )  Rs  denotes  the  serial  resistance,  which  is  a  sum  of  the  membrane 
protonic  and  electronic/ohmic  resistances  and  represents  the 
high-frequency  resistance. 

2)  L  represents  the  high-frequency  cable  inductance. 

3)  The  anode  contribution  contains  Rc t,a  and  Cdi.a,  which  describe 
the  charge-transfer  resistance  of  hydrogen  oxidation  and  the 
capacitance  nature  of  the  anode,  respectively. 

4)  The  cathode  circuit  has  a  resistance  (Rct,c).  which  represents  the 
charge  transfer  of  the  oxygen  reduction.  The  constant  phase 
element  (CPE)  represents  the  double-layer  capacitance  of  the 
real  cathode  with  the  roughness  of  the  catalyst  layer  and  a  non- 
uniform  catalyst  distribution,  and  W  represents  a  finite  length 
Warburg  diffusion  element. 


Rs  L  Ret, a  CPE 


Rs  L  R  „ 

- - qhd-  ■"CD" 

Fig.  1.  Equivalent  electric  circuits  for  Fk/air  (a),  H2/O2  (b)  and  H2/H2  (c)  operating 
conditions. 
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The  EEC  for  the  cathode  electrode  operated  under  air  is  repre¬ 
sented  by  the  well-known  Randles  cell,  which  describes  the  com¬ 
bination  of  kinetic  and  diffusion  processes  in  which  the  capacitance 
element  is  replaced  by  CPE  [83—85]  (Fig.  la). 

The  CPE  represents  the  non-homogeneous  distribution  of 
certain  physical  properties  in  an  electrode,  such  as  the  thickness  or 
roughness  of  the  electrode  surface  that  leads  to  a  distribution  of 
time  constants  in  a  limited  region.  The  impedance  of  the  CPE  was 
introduced  by  Cole  and  Cole  [86]  and  can  be  written  as  [87] 

ZcPE*Qcpe(L)»”  (,) 

where  Qcpe  is  an  admittance  constant  expressed  in  F  cm~2,  pcpe  is 
an  adjustment  parameter  (whose  value  belongs  to  [0, 1],  i  =  (-l)1^ 
2),  and  w  is  an  angular  frequency.  The  pcpe  factor  is  related  to  the 
depression  angle  a  by 

“  =  (1  -  Pcpe)90°  (2) 

If  p  equals  1,  CPE  becomes  a  capacitor: 


Zcpe 


1 

QcpeIw 


(3) 


A  CPE  with  p  =  0.5  can  be  used  to  produce  an  infinite  length 
Warburg  element: 

Zcpe  =  — - 7^=  (4) 

Qcpev  16) 

A  CPE  with  a  p  value  of  0  represents  resistance: 


Zcpe 


1 

Qcpe 


(5) 


The  rate  of  an  electrochemical  reaction  can  be  strongly  affected 
by  the  diffusion  of  reagents  toward  the  catalyst  surface  and  the 
removal  of  products  away  from  the  catalyst.  Whenever  diffusion 
effects  completely  dominate  the  electrochemical  reaction  mecha¬ 
nism,  the  impedance  is  called  the  Warburg  impedance.  Impedance 
of  a  generalized  finite  length  of  a  Warburg  element  can  be  repre¬ 
sented  by  the  following  equation  [35,88]: 


Rw  tanh(i6jTw)Pw 
(iuTw)Pw 


(6) 


The  generalized  finite  length  Warburg  element  is  an  extension  of 
another  more  common  element,  the  finite  length  Warburg.  For  the 
finite  length  Warburg  the  parameter  pw  equals  0.5: 


Zw  = 


Rw  tanhyiwTw 


v/iw7\ 


w 


(7) 


Therefore,  the  finite  length  Warburg  is  characterized  by  resistance 
(Rw)  and  the  time  constant  (Tw): 


Tw 


(8) 


where  <5  is  the  effective  diffusion  thickness  and  D  is  the  effective 
diffusion  coefficient  of  the  particle. 

Using  pure  O2  as  the  oxidant  effectively  eliminates  the  finite 
length  Warburg  diffusion  element  in  a  cathode  circuit  (Fig.  lb)  and 
operation  under  H2/H2  gas  configuration  results  in  the  presence  of 
one  semicircle,  which  might  be  represented  by  EEC  in  Fig.  lc. 


It  was  assumed  that  the  shifts  in  the  impedance  spectra  of 
PEMFC  under  CO  exposure  were  mainly  caused  by  changes  in  the 
anodic  impedance,  which  is  more  complicated  due  to  additional 
reactions  involving  CO  adsorption  and  oxidation  at  the  electro¬ 
catalyst  surface.  The  effects  of  CO  poisoning  can  be  modeled  by 
several  EECs,  as  shown  in  Fig.  2.  The  circuits  (Fig.  2a,  b)  represent 
the  impedance  data  for  electrochemical  reactions  at  the  anode  with 
a  single  adsorbed  intermediate,  where  diffusion  of  the  participating 
species  is  not  rate  limiting. 

In  the  circuit  at  Fig.  2a,  the  capacitance  Ci  is  related  to  the 
steady-state  pseudo-capacitance.  The  resistance,  Ri,  in  parallel  to 
Cj,  is  associated  with  the  adsorbate  formation/removal  37,53,89], 
The  other  resistance,  Rct  a,  is  known  as  the  charge  transfer  resis¬ 
tance  and  relates  to  the  resistance  of  electron  transfer  at  the  anode. 
Together,  these  elements  describe  Faradaic  processes  at  the  elec¬ 
trode.  In  electrochemical  cells,  the  double  layer  capacitor,  Cdi.a,  is 
connected  in  parallel  to  the  Faradaic  impedance  of  the  circuit.  Due 
to  the  dispersion  of  capacitance  over  the  electrode  a  greatly 
improved  fit  of  the  data  might  be  obtained  if  the  double  layer 
capacitor  is  replaced  by  a  CPE.  However,  for  simplification  of  EEC 
model  and  data  fitting  a  pure  capacitance  is  usually  used.  The 
Faradaic  impedance  can  be  expressed  by 


Zf  —  Ret, a  + 


1 

ih 


■  ifoC] 


O) 


The  low-frequency  inductive  loop  is  often  observed  in  imped¬ 
ance  spectra  under  CO  exposure  [36,37,43—46].  The  loop  can  be 
modeled  by  a  pseudo-capacitance  in  parallel  with  a  resistance 
similar  to  Fig.  2a,  and  negative  values  of  the  pseudo-capacitance 
and  resistance  represent  the  inductive  behavior  [37,41,42,53,89], 
A  negative  capacitance  makes  no  physical  sense  in  electrical  terms; 
however,  it  can  be  understood  as  a  decrease  in  the  coverage  of 
adsorbate  with  potential  because 


_  dQ  _  n  4  ^occupied 

~  dE  -  qi  d£ 


(10) 


Fig.  2.  Equivalent  electric  circuits  for  anode  impedance  under  CO  poisoning:  resistance 
in  parallel  to  pseudo-capacitance  circuit  (a),  relaxation  resistance  in  series  with 
inductance  (b),  addition  of  mass  transport  limitations  to  a  previous  case  (c). 
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where  q,-  and  ^occupied  are  the  charge  density  for  the  monolayer  and 
coverage,  respectively.  In  contrast  to  a  negative  C\,  there  is  no  clear 
physical  meaning  of  negative  resistance  Ri.  Negative  resistance 
appears  when  a  current  decreases  with  increasing  voltage;  addi¬ 
tionally,  negative  resistance  is  required  for  a  positive  time  constant 
(r  =  CiRt  )  [41,53],  Harrington  and  Conway  suggested  that  the  fitting 
by  the  inductive  circuit  (Fig.  2b)  is  mathematically  equivalent  to  the 
circuit  in  Fig.  2a  [89],  The  inductive  circuit  is  commonly  used  to 
represent  anode  impedance  for  direct  methanol  fuel  cells  [51—53] 
and  PEMFC  under  CO  poisoning  35,46,49],  The  inductance  LK 
means  that  the  current  signal  follows  a  voltage  perturbation  with  a 
phase  delay  due  to  slowness  of  CO  adsorbed  coverage  relaxation, 
and  Rk  serves  to  modify  the  phase-delay.  In  this  case,  the  Faraday 
impedance  might  be  expressed  by 

Zp=  +  +  (11) 

For  a  general  case,  a  series  combination  of  LK  and  RK  defines  the 
relaxation  impedance,  Zr,  as  it  was  proposed  by  Wagner  et  al. 
[36,38—40],  Conversion  of  CO  is  a  potential-dependent  process 
with  the  rate  k  =  k(r),  which  leads  to 


1  +  itOTK  p  ,  ;  r 

ZK  =  IFd\nk,d£  =  RK  +  t(oLK 

(12) 

Ac  1 

(13) 

K  A/f  JFdln  k/de 

Xr  =  io)TKRK  =  i>’>LK 

(14) 

Here,  I f  denotes  the  Faraday  current,  tr  is  the  time  constant  of 
relaxation,  and  din  fc/dr  is  the  first  derivative  of  the  reciprocal 
relaxation  time  constant  (k  =  1/tr).  The  relaxation  impedance  Zk 
can  be  split  into  the  relaxation  resistance  (Rk)  and  relaxation 
inductivity  (Xr)  with  pseudo-inductance  (Lr  =  trRr)  [38], 

Fig.  2c  presents  an  EEC  where  diffusion  limitations  exist.  The 
circuit  was  suggested  for  a  PEMFC  operated  under  100  ppm  CO 
exposure  and  may  be  expressed  as  [36,38—40] 


Rs  L  Rk  Lk  CPE 


Fig.  3.  EEC  for  EIS  data  obtained  during  operation  under  CO  poisoning  and  various  gas 
configurations:  H2  +  CO/air  (a),  H2  +  C0/02  (b),  H2  +  CO/H2  (c). 


in  a  steady-state  condition  and  can  change  impedance  during 
measurement. 

The  anode  impedance  in  the  Laplace  domain  for  air,  oxygen  or  a 
symmetric  cell  is  described  in  Eq.  (17): 


^anode  (if!l)  — 


Rf(  Ret, a  +  if>LKRci.a 


(to)  Qli.a^ct.af/t  +  !w(Oll,afyc^ct,a  +  Tic)  +  Rk  +  Ret, a 

(17) 


The  cathode  impedance  is  expressed  by  Eq.  (18),  where  the 
impedance  of  the  cathode  double  layer  capacitance  (CPE)  is 
described  by  Eq.  (1 ): 


-  cathode 


(M  =  - 


l 


M'cpeQcpe  +  (V(Rct,c  +  ZWM)) 


(18) 


The  impedance  Zw(i6j)  is  represented  as  a  second-order  function 
[90  .  Eq.  (18)  can  be  simplified  for  pure  O2  operation: 


^cathode  (,0J)  —  . . 


Rr 


(ioj)PcpE  QcpERct,c  + 1 


(19) 


Ret, a  +  Zn 
1  +  Ret, a  /  Zk 


(15)  The  series  impedance  of  the  cell  can  be  calculated  from  the 
following  expression: 


where  Zn  is  the  Nernst  impedance  and  is  equivalent  to  the  finite 
length  Warburg  element  (Zw r).  In  our  work,  the  fuel  cell  operated 
under  2  ppm  of  CO,  which  allows  us  to  assume  that  there  are  no 
diffusion  limitations  of  participating  species  and  that  ZWr  =  0. 
Based  on  this  assumption  and  our  previous  results  [6  ,  the  anode 
EEC  under  CO  poisoning  presented  in  Fig.  2b  was  selected  for  the 
evaluation  and  analysis  of  EIS,  and  the  final  EECs  are  given  in  Fig.  3. 
For  H2/H2  operation  (Fig.  3c),  pure  capacitance  was  transformed  to 
CPE  because  it  provided  more  adequate  fitting  results. 

The  complete  impedance  (Zr)  of  the  cell  under  CO  poisoning  is 
represented  by  the  EECs  shown  in  Fig.  3  and  can  be  grouped  as 
following: 

Zj(oj,  t )  =  Z$((l),  t)  +  2anode(w!  1)  “h  ^cathode(w:  1)  0®) 


Zs(ioj)  =  R$  +  iooL  (20) 


3.2.  Model  for  analysis  of  temporal  evolution  of  current  density 
distribution 

An  attempt  is  made  here  to  apply  a  simple  mathematical  model 
to  achieve  a  qualitative  description  of  the  observed  current  density 
distributions  under  CO  exposure.  During  the  CO  poisoning,  the 
overall  voltage  and  currents  do  not  depend  on  the  frequency 
response  of  the  EECs.  In  a  pseudo-steady  state,  the  impedances 
depend  only  on  resistances  due  to  the  great  values  of  the  time 
constants.  Therefore,  the  complete  impedances  can  be  presented  as 
follows: 


where  Zs  (w,t)  is  the  series  impedance  presented  by  Rs  and  L,  Zan0de 
(oj,t)  is  the  anode  impedance,  and  Zcath0de  («,t)  is  the  cathode 
impedance  contribution.  Moreover,  the  total  impedance  depends 
on  the  applied  frequency  and  time  of  CO  exposure  because  the 
frequency  dependence  is  related  to  the  EEC  dynamics  while  the 
time  dependence  is  associated  with  the  fact  that  the  system  is  not 


Zf(f) 


Rs(t)  +  Rct,c(t)  +  Rw(f)  + 


Rldt)Rct.a(t) 
Rdt)  +  Rct,a(t) 


Z°2(t) 


Rsd)  +  Rct,c(t)  + 


Rl<met,a(t) 

Rl<(t)  +  Ret, a(d 


(21) 

(22) 
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Z"*(t)  =  Rs(t)  +  R(t)  + 


Rl<(t)Rct,a(t) 
Rl<(t)  +  Rct,a(t) 


(23)  =  rT  -  rV 


(37) 


The  dynamic  response  of  the  segment  resistances  (Eqs. 
(21 )— (23))  should  be  represented  by  a  mechanism  that  accounts  for 
the  following  chemical  and  electrochemical  reactions  9,14]: 


H2  +  2[ 

]  <->2[H]  KdH2s 

(24) 

+  H+  +  e-  k+xH 

(25) 

co+  [ ; 

1 — [CO]  Kf°,  /Cs 

(26) 

h2o  + 

[]~[OH]+H+  +  e-  k+f°,  kdlf° 

(27) 

[OH]  +  [CO]  - 2[  }  +  C02  +  H+  +  e-  k+xco 

(28) 

For  this  purpose,  the  CO  poisoning  model  developed  in  Refs. 
[14,91,92]  was  used  here  and  modified  to  describe  the  temporal 
evolution  of  the  spatial  current  densities.  This  model  considers  the 
surface  kinetic  reactions  at  the  anode:  rj  is  the  Tafel  reaction  (Eq. 
(29)),  ry  defines  the  Volmer  contribution  (Eq.  (30)),  raff  is  the  rate  of 
CO  adsorption  (Eq.  (31)),  r^2°  is  the  dissociation  of  H20  (Eq.  (32)), 
and  is  the  CO  oxidation  rate  reaction  (Eq.  (33)). 


Eq.  (36)  accounts  for  the  fact  that  CO  may  be  oxidized  via  a 
chemical  pathway  (rchC0): 

02  +  [  ]  ~  2  [O]  (38) 


[O]  +  [CO]  -  2  [  ]  +  C02 


(39) 


Chemical  oxidation  of  CO  is  controlled  by  the  rate  of  perme¬ 
ation  of  02  through  the  membrane,  and  the  oxidation  kinetic  is 
rapid.  Under  steady-state  diffusion  conditions,  the  flux  of  02  can 
be  expressed  as  [29] 


No2 


I<DPq2 

L 


(40) 


N0l  is  the  steady-state  flux  of  02,  which  crosses  the  membrane 
due  to  diffusion.  K  is  the  Henry's  law  constant,  D  is  the  diffusion 
coefficient  of  02  in  the  membrane,  L  is  the  thickness  of  the  mem¬ 
brane,  and  P02  is  the  partial  oxygen  pressure  in  the  cathode. 
Furthermore,  we  used  the  parameter  p  to  account  for  the  selectivity 
of  02  for  the  oxidation  of  CO  [29]: 


rco 

rch 


4no2 


(41) 


rj  =  K^(p62-K%s82H) 
rv  =  214;%  sinh 


(29) 

(30) 


Considering  the  capacitive  and  faradaic  currents  and  using  the 
equation  of  charge  conservation,  the  time  evolution  of  the  anode 
potential  can  be  calculated  as 

=  €  ~  F(rv  +  r“  +  rdif)  (42) 
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(32) 


rco 


=  kof°dco0OHe 


(33) 


Here,  k„  and  /<„  are  the  reaction  rate  constants  of  the  forward  and 
backward  n  reaction,  Rde2s  and  Kd°s  are  the  equilibrium  constants 
for  H2  and  CO  desorption,  respectively.  The  pressure  in  the 
channel  is  represented  by  p.  T  is  the  average  temperature.  R  and  F 
are  the  gas  constant  and  Faraday  constant,  respectively.  The 
anode  overpotential  is  described  by  rja,  and  the  mole  fraction  of 
CO  is  xco- 

The  balance  of  species  coverage  is  represented  by  Eq.  (34), 
where  the  species  are  CO,  OH  and  H.  Vacant  surface  sites  are  rep¬ 
resented  by  6. 

8  +  0CO  +  #OH  +  #H  =  1  (34) 

The  coverage  dynamics  are  calculated  in  Eqs.  (35)— (37),  and  the 
catalyst  sites  are  described  by  yCt,  where  y  represents  the  rough¬ 
ness  factor  and  Ct  is  the  atom  mole  density  per  cm2. 


FyQ 
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OH  _  H20  _  CO 
(Jf  dis  OX 


(35) 


FyCt 
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co 
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dt 


=  rad  -  or  -  r. 


CO 


ch 


rCO 


(36) 


where  i  is  the  electric  current  density  and  r  is  a  fitting  parameter 
related  to  the  active  area.  The  current  is  the  sum  of  the  faradaic 
current  and  the  capacitive  current.  The  capacitive  current  is  written 
as  the  product  of  the  double-layer  capacity  Cd|  and  the  temporal 
change  of  the  potential  ?ja.  The  faradaic  current  includes  the  H  and 
CO  electrooxidation  currents  (rv  and  )  and  the  net  current  of 
water  dissociation  (r^2°)  [14,92,93], 

To  represent  the  anode  channel  with  a  CO  flow,  a  simple  straight 
channel  has  been  modeled  as  illustrated  in  Fig.  4.  The  model  of  the 
channel  has  been  discretized  into  control  volumes;  each  of  these 
control  volumes  is  connected  to  the  electrochemical  model 
described  above.  The  injected  CO  at  the  channel  inlet  is  dragged  by 
the  anode  gas  flux.  The  effect  of  cross  currents  in  the  GDL  was  not 
taken  into  account  by  the  model,  since  current  density 
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Fig.  4.  Model  diagram  of  straight  channel.  The  elements  TP  and  CV  describe  the 
transport  and  balance  equations,  respectively. 
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distributions,  gradients  of  CO  concentration  and  pressure  drop 
along  flow  channel  are  significantly  larger  than  channel  to  channel. 
The  concentration  of  CO  is  homogeneous  and  constant  in  each 
control  volume  because  well-mixed  gas  conditions  are  assumed. 
The  flux  of  CO  in  the  channel  is  calculated  based  on  the  assumption 
of  pure  advection  of  CO,  and  A,,  depends  on  the  main  flux  in  the 
channel,  in  this  case  the  H2  flux: 

Tco=^coan  (43) 

The  discretization  of  the  channel  allows  us  to  simulate  the  distri¬ 
bution  of  CO  along  the  charnel,  where  the  amount  of  CO  f^yx  to¬ 
ward  the  catalyst  surface  is  j  co.  The  relationship  between  j  co  and 
the  rate  of  CO  adsorption  is  represented  by  the  dimensional 
parameter  8: 

— un  rn 

j  co  =  /5racd°  (44) 

The  total  flux  of  CO  is  the  sum  of  the  CO  transported  in  the  channel 
and  the  consumption  of  CO  in  the  anode  through  the  CO  adsorption 
reaction: 


J  co 


_  dxg?" 

-c^~d ~ 


-UN 

i  co 


(45) 


where  the  capacity  related  to  the  discretized  volume  is  represented 
by  C„. 

The  electric  current  produced  in  each  segment,  ij,  is  calculated 
from  the  following  expression: 

ij  —  (^OC  ~  ?7a ^segment (46) 


10 

'cell  =  E  b  (47) 

J=1 

where  Voc  is  the  open  circuit  voltage  of  the  cell,  and  the  impedance 
of  the  segment  in  steady-state,  ^segment.  is  a  free  variable.  Eq.  (46)  is 
a  simple  polarization  curve  around  a  set  point  ^segment,  where  the 
main  overpotential  is  caused  by  anode  processes.  The  set  point, 
^segment,  is  a  simplified  approach  to  resistance  described  in  Eqs. 
(21)— (23)  depending  on  the  operating  cathode  gas  configuration. 
Eq.  (47)  relates  to  the  total  applied  current,  iceu,  as  a  sum  of  local 
currents,  ij. 

The  model  described  above  of  a  fuel  cell  channel  operated  under 
H2  +  C0/02  conditions  was  simulated  using  the  Modelica  language 
and  Dymola.  For  reducing  the  complexity  of  the  model,  which 
primarily  focuses  on  the  anode  poisoning  processes,  diffusion  and 
flooding  were  not  implemented,  and  the  current  densities  of  all  the 
segments  are  equal  in  the  pre-poisoning  and  recovery  stages.  The 
parameters  used  in  the  simulation  are  listed  in  Table  1  and  were 
acquired  from  the  literature  [14,29,91—94],  with  the  exception  of 
the  values  that  were  not  available,  which  were  fitted  to  reproduce 
the  experimental  data. 

4.  Results  and  discussion 

4.3.  Effect  of  CO  injection  during  H2/O2  operation 

All  the  experiments  were  performed  under  galvanostatic  mode 
with  an  overall  constant  current  density  of  0.8  A  crrT2.  Fig.  5  pre¬ 
sents  profiles  of  the  segment  voltages  and  current  densities 
normalized  to  its  initial  values  vs.  experiment  time.  The  pre¬ 
poisoning  period  was  performed  with  neat  H2/O2  for  70  h.  The 
segment  voltages  and  the  overall  cell  voltage  were  0.665  V  due  to 


Table  1 

Values  of  model  parameters  used  in  the  simulation. 


Parameter 

Description 

Value 

Av 

Advection  gas  flux 

1  1  mirr1 

Qi 

Double  layer  capacitance 

0.450  F 

Ct 

Catalyst  sites  factor 

2.2  x  10~9  mol  cm-2 

Cv 

Gas  capacitance 

2 1 

D 

Diffusion  coefficient  of 

02  in  membrane 

1.7  x  10~6  cm2  s_1 

F 

Faraday  s  constant 

96485.341 5C  mor1 

K 

Henry's  law  constant 

8.86  x  10-7  mol  cm-3  atm-1 

b+tf2 

Kad 

Constant  of  forward  reaction 
rate  of  H2  adsorption 

402  A  cm-2  atm 

i/H2 

Ades 

Equilibrium  constant  for 

H2  desorption 

0.5  atm 

^co 

Constant  of  forward  reaction 
rate  of  CO  adsorption 

50  A  cm-2  atm 

*£ 

Equilibrium  constant  for 

CO  desorption 

1.5  x  10-6  atm 

b+,CO 

Kox 

Constant  of  forward  reaction 
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Fig.  5.  Voltage  (a)  and  normalized  current  densities  (b)  of  individual  segments  vs.  time 
for  an  overall  current  density  of  0.8  A  cm-2  and  2  ppm  CO.  Anode/cathode:  H2/O2, 
0.847/2.018  1  min"1, 100/50%  RH,  48.3/48.3  kPag,  60  °C. 
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the  design  of  the  segmented  cell  system  and  small  size  of  the 
segmented  electrodes  and  cell.  The  performances  of  individual 
segments  were  constant  and  did  not  change  with  time  in  the  pre¬ 
poisoning  phase.  After  the  injection  of  2  ppm  CO  into  the 
hydrogen  stream,  a  fast  decrease  of  the  segments'  voltages 
occurred  and  current  density  redistribution  was  observed.  The 
transition  period  to  reach  steady  state  lasted  approximately  10  h.  At 
steady  state,  the  segments  voltages  settled  at  0.586  V,  which  yiel¬ 
ded  a  performance  loss  of  0.089  V.  Under  CO  exposure,  the  inlet 
segments  1—4  exhibited  the  lowest  current  density  magnitudes, 
while  the  outlet  segments  8—10  showed  an  increase  in  current 
density.  A  decrease  of  the  inlet  segments'  performance  up  to  60%  of 
their  initial  performance  occurred  due  to  adsorption  of  CO  mainly 
at  the  inlet  part  of  the  cell,  which  resulted  in  a  decrease  of  the  CO 
concentration  downstream  and  kept  the  catalyst  layer  effectively 
unpoisoned  at  the  outlet  of  the  cell.  Moreover,  operation  of  the  cell 
in  galvanostatic  mode  required  obtaining  a  constant  current  out  of 
the  entire  cell  and  because  the  current  of  the  inlet  segments 
decreased  after  CO  injection  an  increase  as  high  as  40%  in  the 
current  of  segments  9  and  10  occurred.  A  similar  current  density 
distribution  under  galvanostatic  mode  was  previously  reported  in 
other  papers  [16,22,30];  in  addition,  a  comparable  close  coopera¬ 
tive  behavior  of  Pt  microelectrodes  during  CO  electro-oxidation  was 
presented  in  Ref.  [95],  After  switching  the  fuel  from  the  H2  +  2  ppm 
CO  mixture  to  pure  H2,  the  cell  and  segments  recovered  their  per¬ 
formances  very  close  to  the  initial  voltage  values  of  0.662  V  (vs. 
0.665  V)  within  1  h. 

To  support  the  observed  experimental  trend  in  the  current 
redistribution,  a  simulation  of  the  current  density  was  performed 
using  the  models  described  in  Section  3.2.  Fig.  6a  presents  the 
simulation  results  of  the  spatial  current  density  evolution  with  the 
assumption  that  Eqs.  (24)— (28),  (38),  (39)  represent  the  chemistry 
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Fig.  6.  Simulation  of  current  (a)  and  CO  concentration  (b)  distributions  in  galvanostatic 
mode  under  H2/02  operating  conditions.  The  parameters  used  in  the  simulation  are 
listed  in  Table  1. 


of  the  processes  at  the  anode  and  that  oxidation  of  CO  occurs 
through  electrochemical  and  chemical  pathways.  For  simplification 
purposes,  the  diffusion  and  flooding  processes  were  not  considered 
in  the  model  and  the  current  density  of  all  segments  are  equal  in 
the  pre-poisoning  and  the  recovery  stages.  The  segments  with 
lower  currents  were  located  in  the  inlet,  while  the  outlet  exhibited 
the  greater  performance.  The  dynamic  behavior  of  current  density 
at  the  beginning  of  CO  injection  was  in  good  agreement  with  the 
experimental  data.  The  main  difference  between  the  experimental 
and  simulated  data  is  in  the  dynamics  at  the  recovery  stage.  The 
experiments  showed  a  first  order  decay,  while  the  simulations 
revealed  overshoots  due  to  a  higher  order  dynamics  (Figs.  5  and  6a). 
Time  for  reaching  the  steady  state  at  the  recovery  stage  is  close 
(-4  h)  in  simulations  and  the  real  test.  Moreover,  the  model  predicts 
a  local  CO  concentration  distribution  at  steady-state  (Fig.  6b).  A 
significant  decrease  of  the  CO  concentration  is  observed  for  seg¬ 
ments  7—10.  A  comparison  of  the  simulated  and  experimental  re¬ 
sults  reveals  that  the  proposed  model  adequately  represents  the 
behavior  of  the  fuel  cell  under  CO  exposure  and  might  be  used  for 
prediction  purposes  in  the  future. 

Fig.  7  presents  the  experimental  impedance  data  and  EEC 
modeling  results  for  all  ten  segments  and  the  overall  cell  recorded 
at  different  phases  of  the  experiment.  The  E1S  curve  of  the  fuel  cell 
operated  without  CO  poisoning  consists  of  several  semicircles:  a 
high-frequency  anode  arc  (usually,  it  is  negligible;  however,  it 
might  be  observed  for  low  current  operation),  a  high-frequency 
cathode  arc  attributed  to  a  charge  transfer  resistance  and  double 
layer  capacitance  of  oxygen  reduction  reaction,  and  a  low- 
frequency  arc  presenting  mass  transport  limitations  within  the 
electrode  [32  .  Operation  with  pure  02  as  an  oxidant  eliminated 
mass  transport  obstacles  at  the  cathode  and  simplifies  the  imped¬ 
ance  spectrum,  leaving  a  small  high-frequency  anode  and  signifi¬ 
cant  high-frequency  cathode  arcs,  reflecting  the  charge  transfer 
processes  at  the  electrodes.  These  curves  produced  during  the 
prepoisoning  phase  were  modeled  using  EEC  in  Fig.  lb,  which 
provided  adequate  fitting  to  the  EIS  data. 

After  1  h  of  CO  injection,  there  was  a  slight  increase  in  the 
impedance  of  segments  1—3,  while  the  other  segments  did  not 
show  any  changes.  Moreover,  for  these  segments,  a  drift  of  the 
impedance  values  was  observed  at  low  frequency,  indicating 
changes  in  the  segments'  behavior  during  the  measurement 
because  recording  EIS  at  low  frequency  required  a  finite  time  [54], 
Thus,  we  excluded  the  low-frequency  part  of  the  curves  from  the 
EEC  fitting.  Segments  1—5  exhibited  an  increase  of  the  impedance 
response  after  3  h  of  CO  exposure.  Poisoning  within  5  h  resulted  in 
EIS  changes  for  segments  1—7.  Segment  8  revealed  only  a  slight 
increase  in  EIS  after  49  h,  while  segments  9  and  10  were  not 
affected  by  CO.  Such  a  negligible  impact  of  CO  on  the  outlet  seg¬ 
ments  might  be  explained  by  a  low  CO  concentration  in  the  H2 
steam  (Fig.  6b),  strong  adsorption  of  CO  on  Pt  at  the  inlet  part  of  the 
cell,  and  operation  with  pure  02.  Moreover,  a  typical  pseudo- 
inductive  loop  for  H2  +  CO  operation  was  not  obtained  in  the 
low-frequency  region.  EEC  presented  in  Fig.  3b  was  used  for  fitting 
the  EIS  data  under  CO  exposure.  However,  a  lack  of  pseudo- 
inductive  behavior  allowed  us  to  simplify  the  circuit  to  the  one 
depicted  in  Fig.  lb.  The  EIS  curves  after  1  h  of  recovery  revealed 
significant  decreases  in  the  segments’  impedances.  Total  recovery 
was  observed  after  -4—5  h  of  operation  at  pure  H2. 

High-frequency  region  (10  kHz— 5  kHz)  for  individual  segments 
was  represented  by  pure  inductance  most  likely  due  to  load 
connection  wire  setup  [96],  However,  EIS  curves  for  the  total  cell 
demonstrated  high-frequency  loop,  which  was  close  to  capacitive 
behavior  97—100].  This  behavior  might  be  attributed  to  structural 
properties  of  MEAs,  hardware  size,  measurement  setup  and  cabling 
[96—98],  High-frequency  region  is  not  associated  with  the 
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Fig.  7.  E1S  for  segments  1-10  and  the  overall  cell  during  2  ppm  CO  exposure  test  at  0.8  A  cm-2.  The  experimental  data  are  presented  by  symbols,  and  the  modeling  results  are 
represented  by  solid  lines.  Anode/cathode:  H2/O2, 0.847/2.018 1  min-1, 100/50%  RH,  48.3/48.3  kPag1 60  C.  The  spectrum  at  92  h  was  recorded  up  to  0.01  Hz;  the  others  were  recorded 
up  to  0.05  Hz. 


electrode  processes  and  the  impedance  response  of  the  overall  cell 
exhibited  the  same  trend  as  the  inlet  segments. 

To  understand  the  behavior  of  PEMFC  in  the  presence  of  CO,  the 
mechanism  of  the  CO  and  H2  reactions  should  be  considered.  The 
oxidation  of  hydrogen  on  Pt  was  shown  to  proceed  according  to  the 
Tafel— Volmer  mechanism,  where  the  hydrogen  dissociation  reac¬ 
tion  is  the  rate-controlling  step  [101]: 


2Pt— H  ~  2Pt  +  2H+  +  2e~  (Volmer)  (49) 

Another  mechanism  proposed  for  hydrogen  oxidation  reaction 
(HOR)  is  the  Volmer— Heyrovsky  sequence  [101]: 

H2  +  Pt  «->  Pt-H  +  H+  +  e  (Heyrowsky)  (50) 

Pt— H  <->  Pt  +  H+  +  e~  (Volmer) 


H2  +  2Pt  <->  2Pt— H  (Tafel) 


(48) 


(51) 
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Oxidation  of  hydrogen  on  Pt  occurs  in  the  range  of  0—0.2  V. 
Moreover,  Eqs.  (48)— (51)  indicate  that  the  reaction  is  sensitive  to 
the  surface  coverage  (Pt— H)  and  that  a  decrease  of  available  Pt  sites 
results  in  changes  of  the  reaction  rate.  Moreover,  the  rate¬ 
determining  step  depends  on  the  CO  coverage.  For  dco  <  0.6,  the 
hydrogen  oxidation  reaction  occurs  via  the  Tafel— Volmer  mecha¬ 
nism  with  the  Tafel  reaction  as  the  rate-determining  step  at  the  low 
CO  coverage,  while  the  Volmer  step  controls  the  overall  reaction 
rate  at  the  medium  CO  coverage.  When  the  CO  coverage  is  higher 
than  0.6,  the  Heyrovsky— Volmer  mechanism  is  operative  for  the 
HOR  with  Heyrovsky  sequence  as  the  rate-determining  step  [102], 
Oxidation  of  CO  to  CO2  can  occur  via  chemical  and/or  electro¬ 
chemical  pathways.  The  electrochemical  mechanism  involves  a 
reaction  between  adsorbed  CO  and  oxygen-containing  species,  the 
nature  of  which  remains  under  discussion.  This  species  may  be 
adsorbed  H2O  molecules  or  hydroxyl  (—OH)  species  [103—105],  The 
ignition  potential  of  the  CO  oxidation  on  Pt  at  60  °C  is  in  the  range  of 
0.4— 0.5  V,  and  full  electrochemical  oxidation  occurs  at  0.6— 0.9  V 
[8,27,106,107], 

CO  +  Pt  ^  Pt-CO  (52) 

H20  +  Pt  <->  Pt— H20  (53) 

H20  +  2Pt  ->  Pt-OH  +  Pt-H  (54) 


contributions  of  CO  poisoning,  the  “geometrical  blocking”  or 
“electronic  factor”;  these  factors  should  be  considered  together. 

Chemical  oxidation  of  CO  should  also  be  considered  due  to  the 
diffusion  of  oxygen  from  the  cathode  to  the  anode  compartment.  It 
is  well  known  fact  that  bleeding  a  small  amount  of  oxygen  or  air 
into  the  anode  gas  stream  is  very  effective  in  mitigating  CO 
poisoning  4,13,21—25],  CO  is  oxidized  by  oxygen  via  a  non¬ 
electrochemical  surface  redox  mechanism,  which  is  similar  to  the 
preferential  oxidation  process  [29]: 

02  +  2Pt  <->  2Pt— O  (57) 

Pt-CO  +  Pt-0  -►  2Pt  +  C02  (58) 

The  anode  overpotential  caused  by  the  presence  of  CO  in  this 
case  equals  0.085  V,  which  is  significantly  low  for  the  electro¬ 
chemical  oxidation  of  CO.  Moreover,  no  low-frequency  inductive 
loop  is  present  in  EIS,  which  usually  indicates  a  slow  electro¬ 
chemical  process  involving  a  single  adsorbed  intermediate  forma¬ 
tion/removal  at  an  electrode.  These  reasons  allow  us  to  assume  that 
the  chemical  mechanism  is  dominant,  as  demonstrated  in  our 
previous  paper  [6], 

The  analysis  of  EEC  modeling  of  EIS  data  was  performed  using 
serial  resistance  (R$),  anode  (kct,a)  and  cathode  charge  transfer 
(Rct,c)  resistances.  Fig.  8  shows  the  distributions  of  these  parameters 


CO  H,0  CO  —  H20  —  e~  CO- HO  +  H+  +  e~  CO- 

I  I  I  I  II  I 

Pt  +  Pt  —  Pt  Pt  Pt  Pt  — *  Pt 


CO  HO  CO  -  HO 

II  II 

Pt  +  Pt  Pt  Pt  — ►  CO  2  +  2Pt  +  H+  +  e- 

(56) 


Thus,  CO  more  strongly  chemisorbs  on  the  platinum  sites  than 
hydrogen,  as  indicated  by  the  greater  potential  required  for  the  CO 
oxidation,  and  even  a  small  concentration  of  CO  can  cover  the 
entire  catalyst  surface.  However,  the  rate  of  hydrogen  oxidation  on 
the  few  remaining  platinum  sites  is  rapid  so  it  controls  the  surface 
potential  [8],  Unfortunately,  this  potential  remains  less  than  the 
potential  needed  to  oxidize  CO;  thus,  the  CO  coverage  remains, 
which  is  dictated  by  the  CO  adsorption  isotherm  (Temkin  isotherm) 

[108] ,  The  adsorbed  CO  blocks  sites  for  dissociative  chemisorption 
of  hydrogen,  which  reduces  the  current  because  the  hydrogen 
chemisorption  is  assumed  to  be  second  order  in  the  vacant  sites 
coverage  (8  =  1  -  dco  -  #h  -#oh).  The  coverage  of  CO  is  a  function  of 
the  potential  as  well  as  the  CO  partial  pressure.  In  addition,  CO  and 
H2  do  not  interact  during  their  co-adsorption  at  temperatures  of 
less  than  200  °C.  However,  it  was  noted  that  the  apparent  activation 
energy  for  hydrogen  oxidation  increases  in  the  presence  of  CO 

[109] ,  It  is  very  difficult  experimentally  to  separate  the  two  main 


■HO 

I  (55) 

Pt  — *  C02  +  2Pt  +  H*  +  e- 


as  functions  of  the  segment  positions  at  different  stages  of  the  test. 
Under  CO  exposure,  segments  1—7  exhibited  an  increase  in  the 
serial  resistance,  anode  and  cathode  charge  transfer  resistances.  In 
addition,  segments  8—10  did  not  show  any  changes  in  EIS.  An  in¬ 
crease  in  HFR  for  the  inlet  segments  might  be  due  to  the  local 
current  redistribution,  changes  in  water  production  and  local  water 
management,  which  can  lead  to  insufficient  membrane  hydration. 
As  mentioned  previously,  the  CO  poisoning  of  the  Pt  anode  resulted 
in  an  increase  of  the  apparent  activation  energy  of  hydrogen 
oxidation  and  blocking  of  the  Pt  surface;  thus,  these  geometrical 
and  electronic  effects  were  reflected  by  an  increase  of  the  anode 
charge  transfer  resistance.  The  effective  charge  transfer  resistance 
was  affected  by  the  concentrations  of  the  charged  and  uncharged 
species  (H2O).  The  poisoning  of  the  anode  caused  a  local  inhomo¬ 
geneous  distribution  of  the  generated  H+  and  a  resulting  inhomo¬ 
geneous  distribution  of  the  protonated  and  associated  water 
molecules  at  the  anode  electrode  and  membrane.  This  process 
created  an  uneven  distribution  of  these  species  at  the  cathode, 
which  led  to  an  increase  in  the  cathode  charge  transfer  resistance 
[40].  On  the  other  hand,  CO  poisoning  resulted  in  a  40—60% 
decrease  in  the  current  for  the  segments  1—6  (Fig.  5).  Spatial  EIS 
recorded  for  these  segments  corresponded  to  a  lower  current 
compared  with  the  pre-poisoning  stage.  At  lower  current  densities, 
the  cathode  charge  transfer  resistance  was  larger.  The  EIS  and  EEC 
data  revealed  progressive  CO  poisoning  downstream  as  well  as  a 
dynamic  character  of  the  poisoning  with  the  time  of  the  exposure. 
It  is  clear  that  the  segments  achieved  a  steady-state  after  7—10  h  of 
CO  injection.  The  EIS  curves  after  1  h  of  recovery  showed  significant 
decreases  in  the  segment  impedance  responses,  and  full  recovery 
was  observed  after  -4—5  h  of  stopping  the  CO  injection. 
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Fig.  8.  Distributions  of  serial  resistance  (Rs),  anode  (Rcta)  and  cathode  (Rctc)  charge  transfer  resistances  at  various  CO  exposure  times  and  H2/02  operation.  EEC  fitting  was  per¬ 
formed  in  the  frequency  range  of  10  000-0.05  Hz. 


4.2.  Effect  of  CO  injection  during  H2/H2  operation 

To  study  the  hydrogen  oxidation  reaction  (HOR)  with  CO 
poisoning,  we  performed  the  same  experiment  as  for  the  H2/H2  gas 
configuration.  H2  was  fed  into  the  cathode,  which  served  as  a 
“pseudo”  reference  electrode.  The  hydrogen  for  the  anode  and 
cathode  was  fully  humidified  and  supplied  at  a  constant  flow  rate  of 
0.847  1  min-1. 

The  voltage  responses  and  normalized  current  densities  as  a 
function  of  the  test  time  for  all  the  segments  for  driven  H2/H2 
operation  are  presented  in  Fig.  9.  During  the  pre-poisoning  phase, 
the  segments'  voltages  were  0.130  V.  For  pure  hydrogen  operation, 
the  range  of  individual  segment  current  densities  measured  was 
0.67—0.88  A  cm-2.  After  CO  injection,  the  anode  potential  increased 
to  0.425  V,  resulting  in  a  voltage  difference  of  -0.295  V.  The  current 
density  behavior  was  similar  to  the  previous  case  with  H2/02.  The 
injection  of  2  ppm  CO  resulted  in  a  decrease  in  the  current  densities 
of  the  inlet  segments  (0.51  A  cm-2  for  segment  1)  and  an  increase 
for  the  outlet  segments  (1.16  A  cm-2  for  segment  10).  The  outlet 
segments  had  to  increase  their  currents  to  maintain  the  same 
overall  current  density  of  0.8  A  cm-2.  The  cell  and  segments 
recovered  their  performances  almost  to  their  initial  values  of 
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Fig.  9.  Voltages  (a)  and  current  densities  of  each  segment  normalized  to  the  initial 
performances  (b)  vs.  time  for  an  overall  current  density  of  0.8  A  cnT2  and  a 
contaminant  exposure  of  2  ppm  CO.  Anode/cathode:  H2  +  2  ppm  CO/H2,  0.847/ 
0.847  1  min-1, 100/100%  RH,  48.3/48.3  kPag,  60  °C. 


0.140  V  (vs.  0.130  V)  within  1  h  after  stopping  CO  injection.  The 
observed  difference  in  the  cell  voltage  before  and  after  the  CO  test 
(-10  mV)  might  be  due  to  the  lack  of  oxygen  crossover  and  chemical 
oxidation  of  carbon  monoxide  at  the  anode. 

The  results  of  the  impedance  measurements  and  EEC  fitting  are 
presented  in  Figs.  10  and  11.  It  should  be  noted  that  at  H2/H2 
operation,  EIS  revealed  a  small  semicircle,  which  was  most  likely  a 
combination  of  two  arcs.  The  arc  in  the  high-frequency  region  can 
be  attributed  to  the  process  involving  oxidation  of  the  adsorbed 
hydrogen  atoms  according  to  Eq.  (49),  while  the  arc  in  the  low- 
frequency  region  is  assigned  to  the  dissociative  adsorption  of 
hydrogen  (Eq.  (48))  [31  ].  Alternatively,  the  low-frequency  arc  could 
be  associated  with  the  slow  mass  transfer  processes;  however,  the 
effect  of  a  limited  gas  supply  was  avoided  in  the  work  by  applying 
pure  hydrogen  with  a  high  flow  rate.  As  mentioned  in  Ref.  8],  slow 
dissociation  of  adsorbed  hydrogen  molecules  is  the  rate-limiting 
stage  for  the  oxidation  of  hydrogen.  Thus,  the  EEC  for  fitting  the 
H2/H2  impedance  response  has  only  one  circuit  (Fig.  lc). 

CO  exposure  within  1  h  led  to  an  increase  of  the  impedances  of 
segments  1—3.  Additionally,  for  these  segments,  a  drift  of  the 
impedance  values  was  observed  at  low  frequency,  indicating 
changes  in  the  segments  during  the  measurement.  This  feature  is 
similar  to  that  observed  during  H2/O2  operation.  Injection  of  CO 
within  5  h  resulted  in  a  significant  increase  of  impedance  for  seg¬ 
ments  1  —5  and  a  slight  growth  for  segments  6—9.  The  segments 
reached  steady-state  after  7  h  of  CO  exposure.  The  impendence 
spectra  of  segments  1—8  had  a  negligible  high-frequency  arc  and 
significant  semicircle,  corresponding  to  the  processes  at  the  anode: 
H2  oxidation  in  the  presence  of  CO  and  processes  involving  CO 
conversion.  The  impedance  curves  after  7  h  of  CO  exposure  were 
recorded  up  to  0.05  Hz,  which  allowed  us  to  observe  a  pseudo- 
inductive  loop  for  segments  1  —4.  The  impedance  behavior  under 
these  operating  conditions  can  be  described  by  the  EEC  depicted  in 
Fig.  3c,  where  the  CO  poisoning  affects  oxidation  of  the  adsorbed 
hydrogen  atoms  (Eq.  (49)).  Segments  9  and  10  revealed  more 
complicated  impedance  spectra,  most  likely  due  to  the  effect  of 
water  accumulation  at  the  outlet  part  of  the  cell,  and  they  were 
excluded  from  EEC  fitting.  The  curves  presenting  EIS  after  25  and 
47  h  of  CO  injection  were  measured  in  the  frequency  range  of 
0.01  Hz— 10  kHz.  The  diameters  of  these  curves  were  close  to  those 
of  the  previous  time  point;  however,  the  presence  of  the  pseudo- 
inductive  loop  became  clear  for  segments  1—6.  Segments  7  and  8 
did  not  exhibit  any  inductive  behavior  at  low  frequency.  It  should 
be  noted  that  the  EIS  curves  recovered  after  changing  from  the 
H2  +  2  ppm  CO  mixture  to  pure  H2,  as  did  the  segments'  voltages 
and  current  densities. 
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Fig.  10.  E1S  for  segments  1-10  and  the  total  cell  during  different  CO  (2  ppm)  exposure  phases  at  0.8  A  cnr2.  Anode/cathode:  H2/H2,  0.847/0.847  1  min-1, 100/100%  RH,  60  °C,  48.3/ 
48.3  kPag.  The  experimental  data  are  represented  by  symbols,  and  the  modeling  results  are  represented  by  solid  lines.  The  spectra  at  70  and  92  h  were  recorded  up  to  0.01  Hz;  the 
others  were  recorded  up  to  0.05  Hz. 


Operation  at  H2/H2  gas  configuration  excluded  any  02  pres¬ 
ence,  and  CO  oxidation  should  mainly  go  through  the  electro¬ 
chemical  pathway  (Eqs.  (52)— (56)).  The  anode  overpotential  was 
0.295—0.300  V,  which  was  lower  than  the  ignition  potential  for 
CO  electro-oxidation  on  Pt  at  60  °C.  Yet  even  at  lower  potentials, 
CO  electro-oxidation  is  still  possible  because  the  potential- 
dependent  reaction  rate  constant  of  CO  oxidation  fe(c)  increases 
with  increasing  electrode  potential  31,106];  however,  the  reaction 
rate  might  be  slow.  At  the  beginning  of  CO  exposure  (1—5  h),  the 
anode  potential  was  still  low  for  any  CO  electro-oxidation.  Thus, 


only  CO  adsorption  could  occur  on  Pt,  affecting  oxidation  of  the 
adsorbed  hydrogen  atoms  and  resulting  in  an  increase  of  the 
anode  charge  transfer  resistance  (Rcti  a)  (Fig.  10).  As  soon  as  the 
segment  reached  its  steady-state  and  the  anode  overpotential 
became  -0.300  V,  a  pseudo-inductive  loop  appeared.  The  pseudo- 
inductive  behavior  could  be  attributed  to  slow  CO  electro¬ 
oxidation  and  the  removal  of  CO2  from  the  Pt  surface,  resulting 
in  an  increase  of  the  relaxation  resistance  ( RK )  and  pseudo¬ 
inductance  (Lk)  at  the  beginning  of  the  process  (Fig.  11)  [31,40], 
As  the  electrochemical  oxidation  of  CO  followed  by  desorption  of 


356 


T.V.  Reshetenko  et  al.  /  Journal  of  Power  Sources  2 69  (2014)  344-362 


0.20 


0.12 


30 

1  23456789  10 

25 

0 

CM 

Lw 

■  A 

x  15 

y\ 

3  10 

5 

■M 

0 

123456789  10 

Segments 


-□-44  h  H2 
— 0—46  h,  1  h  CO 
-A-  48  h,  3  h  CO 
-V50h,  5hCO 
52  h,  7  h  CO 
-<—70  h,  25  h  CO 
92  h,  47  h  CO 
-0-165  h,  64  h  H2 


Fig.  11.  Distributions  of  serial  resistance  (Rs),  anode  charge  transfer  resistance  (Rci,.i).  relaxation  resistance  (RK)  and  pseudo-inductance  (Ik)  at  various  CO  exposure  times  for  H2/H2 
operation.  EEC  fitting  was  performed  in  the  frequency  range  of  10  000-0.05  Hz. 


CO2  proceeded,  it  left  the  Pt  surface  vacant  for  hydrogen  chemi¬ 
sorption,  leading  to  a  decrease  of  Rcti  a  and  RK.  This  is  observed  in 
Fig.  11  as  follows  the  initial  increase  of  the  parameters  occurred  at 
the  start  of  poisoning.  Rct,  a,  R/<  and  Lk  began  to  decrease  after 
3—5  h  for  segments  1—4,  while  for  segments  5  and  6  the  decrease 
began  after  7  h. 

4.3.  Effect  of  CO  injection  during  H2/air  operation 

Fig.  12  presents  the  voltage  responses  and  normalized  current 
densities  vs.  experimental  time  for  all  ten  segments  under  H2/air 
operation.  A  pre-poisoning  period  was  conducted  with  pure  H2  for 
the  first  45  h.  The  segment  voltages  and  the  overall  cell  voltage 
were  0.598  V.  It  should  be  noted  that  using  H2/air  and  with  an 
overall  current  density  of  0.8  A  cnrr2,  diffusion  control  of  the 


- Segl, - Seg  2,  Seg  3, - Seg  4,  Seg  5 

- Seg  6,  Seg  7, - Seg  8, - Seg  9, - Seg  10 


Fig.  12.  Voltage  (a)  and  normalized  current  densities  (b)  of  individual  segments  vs. 
time  for  an  overall  current  density  of  0.8  A  cnT2  and  2  ppm  CO.  Anode/cathode:  H2/air, 
0.847/2.018  1  min-1, 100/50%  RH,  48.3/48.3  kPag,  60  °C. 


reaction  in  the  downstream  segments  is  more  dominant  vs  H2/02 
and  H2/H2  operation,  especially  for  the  outlet  segments  9  and  10, 
which  suffered  from  water  flooding  as  well  as  the  low  02  partial 
pressure  in  the  gas  stream.  Such  performance  distribution  at  the 
pre-poisoning  step  was  confirmed  by  spatial  EIS.  The  injection  of 
2  ppm  CO  into  the  H2  stream  resulted  in  a  fast  decrease  in  the 
segments'  voltages  and  current  redistribution  (Fig.  12).  After 
7—10  h,  the  cell  reached  a  steady-state  with  a  voltage  of  0.320  V. 
The  voltage  drop  was  -0.280  V,  and  the  current  redistribution 
varied  in  the  range  from  -40%  to  32%  from  the  initial  segments' 
performance.  Such  changes  are  attributed  to  progressive  CO 
poisoning  of  the  catalyst  at  the  inlet  segments  during  the  first  hours 
of  CO  injection.  The  catalyst  poisoning  was  also  accompanied  by  a 
reduction  of  the  electrochemical  reaction  rate  and  the  current 
produced;  thus,  the  current  density  of  the  inlet  segments  decreased 
significantly.  Because  the  cell  was  operated  under  galvanostatic 
control  and  the  overall  cell  current  was  maintained  constant,  there 
was  a  current  increase  for  the  outlet  segments  to  compensate  for 
the  losses  at  the  inlet  segments.  After  changing  the  anode  gas  from 
the  H2  +  2  ppm  CO  mixture  to  pure  H2,  the  cell  and  segments' 
voltages  recovered  to  0.595  V,  and  the  segments'  currents  became 
similar  to  their  initial  values. 

Fig.  13  presents  the  impedance  spectra  for  all  ten  segments,  the 
overall  cell  during  different  stages  of  the  experiment,  and  the 
results  of  the  EEC  fitting.  A  comparison  of  the  curves  recorded 
under  neat  H2  operation  revealed  the  existence  of  a  small  high- 
frequency  anode  loop,  a  cathode  loop,  and  a  low-frequency  arc 
corresponding  to  02  mass  transfer  limitations  at  the  cathode. 
There  was  an  increase  in  the  diameter  of  the  low-frequency  loop 
from  the  inlet  to  the  outlet,  reflecting  a  growing  diffusion  limi¬ 
tation  downstream. 

A  comparison  of  EIS  curves  under  CO  exposure  revealed  an  in¬ 
crease  in  impedance  for  segments  1—3  after  1  h,  while  the  others 
did  not  exhibit  any  significant  changes.  Such  changes  were  attrib¬ 
uted  to  progressive  CO  poisoning  of  the  catalyst  at  the  inlet  seg¬ 
ments  during  the  first  hours  of  the  CO  injection.  Additionally,  for 
these  segments,  a  drift  of  the  impedance  values  was  observed  at 
low  frequency  during  recording  of  the  EIS  curves;  this  feature  was 
similar  to  that  observed  during  H2/02  and  H2/H2  operation. 
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Fig.  13.  EIS  for  segments  1-10  and  the  total  cell  during  different  CO  (2  ppm)  exposure  phases  at  0.8  A  cm-2.  Anode/cathode:  H2/air,  0.847/2.018  1  min-1, 100/50%  RH,  60  °C,  48.3/ 
48.3  kPag.  The  experimental  data  are  represented  by  symbols,  and  the  modeling  results  are  represented  by  solid  lines.  The  spectra  at  115  and  140  h  were  recorded  up  to  0.01  Hz,  and 
the  others  were  recorded  to  0.05  Hz. 


Based  on  the  EIS  data,  segments  1  —6  reached  steady-state  after 
5  h  of  CO  injection.  There  was  an  increase  in  HFR,  the  anode  and 
cathode  loops  for  the  inlet  segments.  Additionally,  segments  1—6 
had  pseudo-inductive  loops  in  the  low-frequency  region,  which 
reflects  the  existence  of  the  slow  process  of  CO  oxidation  involving 
at  least  one  intermediate.  Segments  7—9  reached  steady-state 
several  hours  later.  Segment  10  exhibited  an  increasing  low- 
frequency  loop,  which  finally  became  negative  and  was  located  in 
quadrant  II  of  the  Nyquist  plot  (Fig.  13).  Such  EIS  behavior  made  EEC 


fitting  difficult,  and  segment  10  was  excluded  from  analysis.  The 
other  segments'  EIS  curves  could  be  described  by  the  EEC  depicted 
in  Fig.  3a.  EIS  after  70  and  95  h  of  CO  exposure  were  recorded  in  the 
frequency  range  of  0.01  Hz— 10  kHz.  The  presence  of  the  pseudo- 
inductive  loop  became  pronounced  for  segments  1—6.  Segments 
7—9  did  not  demonstrate  any  inductive  behavior  at  low  frequency. 

Experimental  EIS  measurements  were  performed  with  a  time 
interval  of  several  hours  (2—5  h);  thus,  it  was  beneficial  to  evaluate 
the  EIS  evolution  in  detail  for  the  first  10  h  of  CO  exposure.  For  this 
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purpose,  it  was  necessary  to  reduce  the  time  interval  between  the 
measurements;  however,  it  was  impossible  to  perform  EIS 
recording  so  frequently,  since  a  certain  amount  of  time  was 
required  (-0.5—1  h).  To  resolve  this  obstacle,  cubic  interpolation  of 


the  experimental  results  was  applied.  EIS  interpolation  was  per¬ 
formed  in  the  time  range  from  0  to  40  000  s  (0—11.11  h)  with  an 
interval  of  1000  s  (0.28  h).  The  EIS  data  of  segment  10  are  not 
presented  here  because  the  interpolation  was  not  performed  due  to 


Fig.  14.  Time  evolution  of  EIS  for  segments  1-9  using  cubic  interpolation  performed  in  the  frequency  range  of  10  000-0.05  Hz. 
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the  complicated  nature  of  the  recorded  spectra  (Fig.  13).  The 
interpolated  data  profiles  are  presented  in  Fig.  14  and  show  the 
detailed  temporal  evolution  of  the  impedance  behavior,  which 
could  not  be  detected  by  the  experimental  measurements.  A 
gradual  increase  of  the  total  impedance  of  the  segments  and  the 
development  of  the  low-frequency  pseudo-inductance  can  be 
clearly  distinguished.  A  combination  of  the  proposed  temporal 
current  density  evolution  model  with  the  EIS  interpolation 
approach  can  be  a  powerful  tool  for  detailed  analysis  of  fuel  cell 
performance.  The  method  can  be  applied  for  studying  the  effects  of 
different  operating  conditions  and  fuel/air  contaminants  and  for 
the  prediction  of  fuel  cell  performance. 

EEC  fitting  revealed  that  under  CO  exposure,  there  were  in¬ 
creases  in  HFR,  anode  and  cathode  charge  transfer  resistances,  as 
well  as  in  the  cathode  Warburg  resistance  (Fig.  15).  Analysis  of  the 
segments’  EIS  parameters  distributions  showed  that  the  anode  and 
cathode  charge  transfer  resistances  decreased  downstream,  and 
the  effect  of  CO  on  the  outlet  segments  decreased  due  to  the  con¬ 
sumption  of  CO  by  previous  segments  and  the  reduction  of  CO 
concentrations  downstream.  The  effect  of  CO  on  segment  10  was, 
therefore,  most  likely  negligibly  small,  and  its  EIS  behavior  may  be 
explained  by  mass  transfer  limitations  and  02  depletion  [63,64] 
instead  of  CO  poisoning. 

Operation  at  F^/air  configuration  most  likely  resulted  in  CO 
oxidation  through  chemical  and  electrochemical  pathways  (Eqs. 
(52)— (58)).  The  anode  overpotential  was  0.280—0.285  V,  which 
was  lower  than  that  for  H2/H2  and  the  ignition  potential  for  CO 
electro-oxidation  on  Pt  at  60  °C.  However,  even  at  lower  potential, 
slow  CO  electro-oxidation  is  possible  as  demonstrated  previously 
[31,106].  Additionally,  oxygen  can  diffuse  from  the  cathode  to  the 
anode  and  contribute  to  the  chemical  CO  oxidation  (Eqs.  (57)  and 
(58))  as  discussed  earlier.  At  the  beginning  of  CO  exposure 
(1—5  h),  adsorption  of  CO  occurred  on  Pt,  affecting  the  oxidation  of 
the  adsorbed  hydrogen  atoms  and  generated  current  density, 
which  caused  an  increase  of  the  anode  charge  transfer  resistance 


(Rct,  a)  and  HFR  (Fig.  15)  [40],  Carbon  monoxide  also  affected  the 
cathode  and  led  to  an  increase  in  the  cathode  charge  transfer 
resistance  (Rct,  c)  due  to  uneven  distributions  of  protonated  and 
associated  water  molecules  at  the  anode,  membrane  and  the 
cathode  [12,40],  In  addition,  there  was  an  increase  of  mass  trans¬ 
port  resistance  (Rw).  After  5  h  of  CO  exposure,  a  pseudo-inductive 
loop  was  observed,  which  could  be  connected  to  slow  CO  electro¬ 
oxidation  and  removal  of  C02  from  the  Pt  surface,  resulting  in  an 
increase  of  relaxation  resistance  (Rk)  and  pseudo-inductance  (Lk) 
(Figs.  13  and  15).  As  the  electrochemical  oxidation  of  CO  proceeded, 
a  reduction  in  Rk  at  a  steady-state  was  also  observed  (Fig.  15) 
[31,40],  A  comparison  of  Rcti  a  behavior  under  the  H2/H2,  H2/02  and 
H2/air  operating  conditions  (Figs.  11  and  15)  revealed  that  the  value 
of  Rct,  a  at  H2/air  is  lower  than  that  at  H2/H2  and  higher  than  that  of 
H2/O2  which  could  be  attributed  to  a  contribution  of  the  chemical 
CO  oxidation  to  the  carbon  monoxide  removal  from  the  Pt  anode 
surface.  EIS  of  the  total  cell  exhibited  the  same  trend  that  was 
observed  for  the  individual  inlet  segments:  an  increase  in  the 
impedance  response  and  existence  of  a  pseudo-inductive  loop  at 
low  frequency. 

After  changing  the  anode  gas  from  the  H2  +  2  ppm  CO  mixture 
to  pure  H2,  the  cell  segments’  EIS  curves  demonstrated  significant 
decreases  in  the  impedance  responses,  and  the  curves  returned  to 
their  initial  states  after  several  hours  of  recovery. 

4.4.  Effect  of  CO  exposure  on  fuel  cell  performance  and  ECA 

Fig.  16a— c  presents  the  voltage  distributions  at  fixed  current 
densities  as  functions  of  the  segment  locations  before  and  after 
the  CO  exposure  tests  performed  for  the  different  gas  configu¬ 
rations:  H2/air,  H2/O2  and  driven  H2/H2.  A  comparison  of  the 
curves  revealed  that  after  the  first  CO  exposure  test  under  H2/air 
conditions,  the  performance  did  not  change  significantly 
(Fig.  16a).  There  was  a  slight  decrease  in  the  performances 
(0.015—0.030  V)  for  segment  2  mainly  due  to  ohmic  losses. 


Segments  Segments 


20  h  H2;-0— 46  h,  1  h  CO; -A- 48  h,  3  h  CO;  -O“50  h,  5  h  CO 

115  h,  70  h  CO;  140  h,  95  h  CO;  146  h,  1  h  H2;  164  h,  19  h  H2 


Fig.  15.  Distributions  of  serial  resistance  ( Rs ),  anode  (Rct,a)  and  cathode  (Rc t,c)  charge  transfer  resistances,  Warburg  resistance  (Rw),  relaxation  resistance  ( RK )  and  pseudo-inductance 
( Lk )  at  various  CO  exposure  times  for  H2/air  operation.  EEC  fitting  was  performed  in  the  frequency  range  of  10  000-0.05  Hz. 
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Fig.  16.  Distributions  of  electrochemical  properties.  Cell  voltage  distribution  for  various  current  densities  up  to  1.2  A  cur2  before  and  after  exposure  to  2  ppm  CO  and  different 
oxidant  gases:  air  (a),  02  (b)  and  H2  (c).  Anode/cathode  for  performance  measurements:  H2/air,  2/2  stoichiometry,  100/50%  RH,  60  °C,  48.3/48.3  kPag.  Anode  and  cathode  ECA 
differences  between  after  and  before  CO  exposure  tests  performed  with  air  (d),  02  (e)  and  H2  (f)  as  oxidants.  The  anode/cathode  parameters  for  ECA  measurements:  H2/N2,  0.75/ 
0.75  1  min-1, 100/100%  RH,  35  °C,  0/0  kPag,  20  mV  s-1. 


Fig.  16b  shows  the  voltage  distributions  before  and  after  the 
second  CO  exposure  test  under  H2/O2  conditions.  The  curves 
indicated  that  after  the  second  CO  test,  there  was  a  slight  per¬ 
formance  decrease  (0.015—0.020  V)  in  inlet  segments  1—4,  while 
the  others  had  the  same  performances.  An  analysis  of  the  spatial 
performance  before  and  after  the  third  CO  test  (H2  +  2  ppm  CO/ 
H2)  revealed  a  slight  increase  in  the  performances  of  segments 
1—5  (-0.025  V)  over  the  entire  current  range  (Fig.  16c).  The 
observed  slight  performance  variations  demonstrated  that  CO 
exposure  did  not  damage  the  MEA. 

The  distributions  of  the  anode  and  cathode  ECA  differences  are 
shown  in  Fig.  16d— f.  The  evolution  of  ECA  after  the  2  ppm  CO  +  H2/ 
air  test  indicated  that  the  anode  and  cathode  losses  were  6.40  and 
7.11%,  respectively.  It  should  be  noted  that  such  an  ECA  decrease 
was  observed  during  fuel  cell  operation  without  any  CO  exposure 
and  might  be  related  to  typical  degradation  processes  in  occurring 
in  new  MEAs.  After  the  CO  poisoning  under  H2/02  operating  con¬ 
ditions,  a  decrease  of  the  overall  cathode  ECA  was  13.53%,  while  the 
anode  ECA  decrease  was  only  2.49%.  A  similar  ECA  trend  was  pre¬ 
viously  observed  in  our  paper  [6],  The  CO  exposure  under  H2/H2 
driven  test  conditions  resulted  in  a  slight  increase  of  anode  and 
cathode  ECA,  which  allowed  us  to  conclude  that  ECA  was  not 
seriously  affected  during  the  test.  After  the  CO  test,  there  was  no 
pinhole  formation  in  the  MEA,  and  the  overall  hydrogen  crossover 
current  was  -0.63  mA  citT2. 

5.  Conclusions 

The  effects  of  low  CO  concentration  (2  ppm)  close  to  the  H2  fuel 
specification  (0.2  ppm,  ISO  14687-2,  2012)  on  spatial  PEMFC  per¬ 
formance  at  a  constant  overall  cell  current  under  H2/air,  H2/O2,  H2/ 
H2  conditions  were  studied  using  a  segmented  cell  system  and 
spatial  EIS.  CO  injection  into  H2  stream  caused  a  significant  current 
redistribution.  The  current  densities  of  inlet  segments  1—4 
decreased,  while  outlet  segments  7—10  demonstrated  an  increase 


in  currents.  In  addition,  the  cell  voltage  decreased  depending  on  the 
applied  cathode  gas  (air,  02,  or  H2).  A  mathematical  model  from 
literature  [9,14,29,92]  was  adapted  to  describe  and  understand  the 
temporal  evolution  of  the  current  density  distribution  under  CO 
poisoning.  The  proposed  simple  model  exhibited  an  adequate  dy¬ 
namic  response  of  local  currents  during  CO  poisoning.  This  suc¬ 
cessful  revelation  makes  the  further  refinement  and  development 
of  the  model  very  promising  for  its  potential  application  to  study 
the  effects  of  other  fuel/air  impurities  on  the  spatial  fuel  cell 
performance. 

Spatial  EIS  data  profiles  were  analyzed  by  EEC  modeling.  The 
results  showed  an  increase  in  the  HFR,  anode  and  cathode  charge 
transfer  resistances  under  CO  exposure.  Operation  with  O2  as  an 
oxidant  did  not  result  in  the  appearance  of  a  pseudo-inductive  loop 
at  low  frequency;  moreover,  the  anode  overpotential  was  only 
0.085  V,  which  is  lower  than  the  ignition  potential  for  CO  electro¬ 
oxidation  on  Pt  at  60  °C.  These  observations  allowed  us  to  as¬ 
sume  that  the  oxidation  and  removal  of  adsorbed  CO  occurred 
predominantly  through  a  chemical  pathway,  where  diffused  O2  can 
directly  react  with  CO.  In  contrast,  operation  under  the  H2/H2 
driven  mode  led  to  higher  anode  overpotential  (-0.300  V)  and  the 
existence  of  the  pseudo-inductive  loop  for  segments  1—6.  The 
pseudo-inductive  behavior  might  be  attributed  to  a  surface  relax¬ 
ation  process  of  the  anode  including  slow  electrochemical  reactions 
(CO  electro-oxidation)  with  a  single  adsorbed  intermediate.  The 
lack  of  oxygen  led  to  CO  oxidation  and  removal  mainly  via  an 
electro-oxidation  mechanism.  Using  air  as  the  feed  gas  at  the 
cathode  provided  conditions  for  CO  oxidation  through  the  combi¬ 
nation  of  chemical  and  electrochemical  pathways,  since  the 
pseudo-inductance  for  segments  1—6  was  observed,  and  the  anode 
overpotential  was  high  enough  for  slow  CO  electro-oxidation 
(-0.280  V). 

A  comparison  of  the  EEC  parameters’  distributions  revealed  a 
decrease  of  Rctia,  Rct,c.  i/c  with  the  segments’  location.  This  trend 
in  the  parameters  distributions  together  with  the  lack  of  any 
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pseudo-inductive  behavior  for  the  outlet  segments  7—10  under  H 2/ 
air  and  H2/H2  operating  conditions  allowed  us  to  assume  that  CO 
adsorption  and  Pt  poisoning  mainly  occurred  at  the  inlet  part  of  the 
cell  within  5—10  h,  providing  dynamic  responses  of  the  segments 
and  their  spatial  performance.  A  combination  of  the  current  dis¬ 
tribution  model  with  an  EIS  interpolation  approach  was  demon¬ 
strated  to  be  a  powerful  tool  for  detailed  analysis  and  prediction  of 
fuel  cell  spatial  performance. 
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